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00 . Abstract 

O 

We calculate the vacuum polarization functions on the lattice using the overlap fermion formu- 
/\ ' lation. By matching the lattice data at large momentum scales with the perturbative expansion 

supplemented by Operator Product Expansion (OPE), we extract the strong coupling constant 
as{ii) in two- flavor QCD as ^\fq — 0.234(9) (_ q) GeV, where the errors are statistical and sys- 
tematic, respectively. In addition, from the analysis of the difference between the vector and 
axial-vector channels, we obtain some of the four-quark condensates. 
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I. INTRODUCTION 

In Quantum Chromodynamics (QCD) the vacuum polarization, defined through the 
(axial-)vector current correlator, contains rich information of its perturbative and non- 
perturbative dynamics. In the long distance regime it is sensitive to the low-lying particle 
spectrum. The short distance regime, on the other hand, can be analyzed using perturbation 
theory supplemented by the Operator Product Expansion (OPE). The current correlator can 
be expressed as an expansion in terms of the strong coupling constant as together with power 
corrections of the form {O^^'^) / Q^ . Here, the local operator O^^'^ has a mass dimension n and 
Q is the momentum scale flowing into the correlator. Determination of a^ (and of the vac- 
uum expectation values {O^^'^), in principle) can be performed by applying the formulae for 
experimental results of e~^e~ cross section or r decay distributions l|, for instance. On the 
other hand, if one can calculate the correlators non-perturbatively, theoretical determination 
of those fundamental parameters is made possible. 

Lattice QCD calculation offers such a non-perturbative technique. Two-point correlators 
can be calculated for space-like separations. In this work we investigate the use of the 
perturbative formulae of the correlators for the lattice data obtained in the high Q^ regime. 
The strong coupling constant ag may then be extracted. In such an analysis, it is essential 
to find the region of Q^ where the perturbative expression can be applied and at the same 
time the discretization error is under control. By inspecting the numerical data, we find 
that this is indeed possible at a lattice spacing a ^ 0.12 fm if we subtract the bulk of the 
discretization effects non-perturbatively. The remaining effect can be estimated using the 
perturbation theory. 

The idea of analyzing the short distance regime is not new: in fact, the analysis of hadron 
correlators in the whole length-scales was proposed 15 year ago |3], but to our knowledge 
quantitative analysis including the determination of a^ and (O*^")) has been missing until 
recently. (Calculation of the vacuum polarization from the vector current correlator in lattice 
QCD may be found in 3|, |j|. More recently, an analysis of charmonium correlator has been 
pubhshed |5|.) 

While the vacuum polarizations Ilj{Q^) (J denotes vector or axial-vector channel) are 
ultraviolet divergent and their precise value depends on the renormalization scheme, their 
derivative Dj{Q^) = —Q'^dIlj{Q'^)/dQ'^, called the Adler function 6(], is finite and renormal- 



ization scheme independent. Therefore, the continuum perturbative expansion of Dj{Q^) to 
order a^ [7|, |8| , can be directly apphed to the lattice data. At relatively low Q'^ region, higher 
order terms of OPE become relevant. They include the parameters describing the gluon con- 
densate {agG"^) and the quark condensate {mqq) (we suppress quark flavor index assuming 
degenerate up and down quark masses) at 0{1/Q^), and four-quark condensates (Og) and 
(Oi) at 0(l/g^) y, \l^. (The exphcit form of Og and d will be given in Section [11111) 
We use the lattice QCD data containing two dynamical flavors described by the overlap 



fermions 



11| . The simulations are performed at lattice spacing a = 0.118(2) fm on a 16^ x 32 



lattice. For the details of the simulation including the choice of the lattice actions and 



parameters, we refer llj. The physical volume is about (1.9 fm)^, which is relatively small 
compared to the present large scale QCD simulations. The finite volume effect is, however, 
not significant for the short distance quantities considered in this work. The quark masses 
mq in this analysis are 0.015, 0.025, 0.035 and 0.050 in the lattice unit, that cover the range 
[ms/Q,ms/2] with rris the physical strange quark mass. An analysis of pion mass and decay 



constant is presented in 



V. 



The main advantage of this data set is that both the sea and valence quarks preserve exact 



chiral and flavor symmetries by the use of the overlap fermion formulation 13|, [ij] . (Although 
the fermionic currents used in our calculation are not conserved at finite lattice spacings, it 
does not change the following argument of the operator mixing.) The perturbative formulae 
for the vacuum polarizations can therefore be applied without any modification due to 
explicit violation of the chiral symmetry. For instance, the scalar density operator qq to 
define the quark condensate is free from the leading power divergence which scales as 1/a^. 
This means that a term of the form ma~'^ /Q^ is forbidden in the OPE formula as in the 
continuum theory. With the Wilson-type fermion formulation, this term may appear and has 
to be identified and subtracted non-perturbatively. With the staggered fermion formulation, 
there is no such problem because of its remnant chiral symmetry, while the effect of taste- 
breaking may become significant when {aQY becomes 0(1). 

This paper is organized as follows. In Section [TTl we define the vacuum polarization 
functions and explain the method to calculate them on the lattice. Subtraction of lattice 
artifacts is discussed in some detail. Section IIIII summarizes the perturbative formulae of 
OPE. Then, in Section [IV] we show the results of fitting of our data with the perturbative 
formulae. Estimate of the systematic errors is also given. Conclusions are given in Section IVl 



II. VACUUM POLARIZATION FUNCTION 



Definition 



In the continuum theory, the vacuum polarization functions IIj (Q^) are defined through 
two-point correlation functions as 



{J,J:){Q) = j d'xe^Q-{T{r^{x)Ji\Q)}) 

where the current J*-' may either be a vector current V*-' = (ji'jfj.qj or an axial-vector current 
^u = Qi'liJ,'l5Qj with fiavor indices i ^ j. IIj (g^) and IIj (g^) denote the transverse and 
longitudinal parts of the vacuum polarization, respectively. For the vector channel {J = V), 
Ily (g^) = is satisfied due to current conservation. For the axial- vector channel {J = A), 
the longitudinal component may appear when the quark mass is finite. 

In the lattice calculation we employ the overlap fermion formulation [13|, [iJ] , for which 
the Dirac operator is given by 

D{m) = {niQ + yj + (jrio - yj TsSgn [Hw{-mo)] (2) 

for a bare quark mass m. The kernel operator Hw{—mo) = 75-Dh/(— mo) is constructed from 
the conventional Wilson-Dirac operator D]^{—mQ) at a large negative mass —rriQ. We set 
rriQ = 1.6 in the numerical simulation. We use the vector and axial-current operators of the 
form 



D 

2mo 



V« = Zqa, ( 1 - ;^ ) %. (3) 



^ = ^^^-r. (i - ^) *. (4) 

With this choice, the vector and axial charges form a multiplet under the axial transformation 
S'Xqi = er^j'y^^l — D/mQ)qj, 5\qi = sqjT'^f^^, where e denotes an infinitesimal parameter and 
t"- is a generator of the fiavor SU(2) symmetry. The overlap fermion action is invariant 
under this modified chiral transformation lisl, as it satisfies the Ginsparg- Wilson relation 

n 

-D75 + 75-D = Dj^D/niQ [16|. The common renormalization factor Z has been calculated 
non-perturbatively as Z = 1.3842(3) [l 



An obvious drawback of the (axial-)vector currents in ([3]) and (j3]) is that the current 
conservation property dfj,J^ = {J = V or A) is not satisfied at finite lattice spacing. It 
leads to a significant complication in the extraction of the functions Ilj (Q^) and II j (Q^) 



as described in the next subsection. The use of the conserved (axial-) vector current 17 1 
reduces this complication. Once we have extracted the functions Ily (Q^) and IIj (Q^), 
these two types of currents should give an equally good approximation to the continuum 
one up to the unphysical constant shift (and the discretization error). Our preliminary study 
employing the conserved currents shows that this is indeed the case. 

B. Non-perturbative subtraction of lattice artifact 

Due to the discretization effects including the current non-conservation effect, the two- 
point correlation functions ([I]) may have more complicated structures. Taking account of 
remaining symmetries on the lattice (parity and cubic symmetries) but without the current 
conservation, the correlators on the lattice (JfiJuY^^iQ) can be expressed as an expansion in 

Q,- 

oo oo 

- J2^niQ)QTs,u - Y. cLiQ){Ql"'^'Q'r' + Q'r-^'Qf-'}, (5) 

n=0 m,n=l 

in the momentum space. The lattice momentum Q^ is defined as Q^ = (2/a) sin(7m^/L^) 
with an integer four-vector n^ whose components take values in (— L^/2, L^/2] on a lattice 
of size Lfj^ in the fi-th direction (Lj=i^2,3 = 16 and Lt = 32 in our case). The functions 
corresponding to the continuum counterparts, Ilj (Q) and 11 j (Q) (= Hj (Q) -t-Ilj (Q)), 
may also have Lorentz-violating effects and could be a function of Q^ in general rather than 
a function of just a single argument Q^. 

The term Bq{Q)5^jj^ which has the same Lorentz structure as the term of physical lij {Q) 
does, contains a quadratically divergent contact term. Since one cannot disentangle the 
physical contribution from the unphysical divergence using the Lorentz structure alone, we 
focus on extracting 11 j (Q), which is free from the contact term. 

The terms including functions Bi^^Q^Q) and C;^„((5) represent the lattice artifacts that 
violate the Lorentz symmetry. They are generally written in terms of an expansion in aQ^ 
and aQ^. (Physically relevant terms are separately written with a conventional notation 



Uj (Q) and njy (Q)-) The lowest order term B({Q) remains constant in the continuum 
hmit aQ -^ 0, while the terms of B2{Q) and C(^{Q) are relatively suppressed by 0{{aQY) 
and vanish in the continuum limit. Higher order terms are suppressed by additional powers 
of a at a fixed Q. Since the momentum scale Q of interest is not much less than the 
lattice cutoff 1/a, the convergence of the expansion at our lattice spacing must be carefully 
investigated for the lattice data. These terms can be identified non-perturbatively, and we 
found that the lowest non-trivial terms including i?^((5) and C(i{Q) are already very small 
as described below. Higher order terms are thus safely neglected. 

Extraction of -B/glQ) ^"^^ ^iiiQ) from the lattice data goes as follows. The off-diagonal 
components {J^iJuY^^iQ) (/i 7^ i^) contain Hj (Q) and C/^((5), hence by taking the data 
with two different momentum configurations giving the same Q^ one can solve a linear 
equation to disentangle Hj (Q) from the lattice artifact. To be explicit, for two different 
momentum configurations aQ^^^ and aQ^'^^ giving the same {aQ^^^Y = {aQ^'^^Y = {clQYi the 
linear equation is written as 

{JM'XMQ^'^) = aQ?^«Qi^^nr^)(g(^))- (ag(^)(ag(^))^ + ag(^)(ag(^))^)c/,(g(^)). 

(6) 

We may assume the equalities U^j^^\q^^^) = U^j^^\q^^^) and C^^{Q^^^) = C(^{Q^^^) for 
small enough (aQ)"^, because aQ^^^ and aQ^"^^ are different only by permutations of space- 
time directions. The linear equation ([6]) can be solved when 

Qi^Q^^^Q?^Q^^^ [iQ?^r + {Q^J'r - {QPr - (Q^J^n + 0. (7) 

It is easy to see that three different non-zero components must be contained in aQ^^^ and 
aQ^"^^ to satisfy ([7]). The smallest possible momentum assignment corresponds to the combi- 
nation |n/j I = (2,1,0,1), |n/i I = (1,2,0,1) with (/i, z^) = (1,4) and its permutations. Since 
the fourth (temporal) direction is longer for our lattice (i^j=i,2,3 = 16 while L4 = 32), the 
fourth component of Q^ is effectively 1/2 of spatial components when they are the same 
in n^. The corresponding momentum squared for this choice is [aQ^ ~ 0.776. For larger 
lattice momenta, there are many possible choices that this procedure is applied. 

The lattice artifact in the diagonal pieces, B({Q) and B2{Q), can be extracted in a 
similar manner by solving linear equations for n = u after subtracting the C^i{Q) terms. 
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FIG. 1: Momentum dependence of B({Q), B({Q)/o?, and C(^{Q)/a^ at 



nin 



0.015. Circles 



(crosses) show the vector (axial-vector) channel. The solid curves represent a polynomial fit and 
the dashed curves show the one-loop results. 



For instance, the leading contribution {lij {Q)Q'^ — Bq{Q))5^u is extracted by subtracting 
the sub-leading contribution B({Q)Q'i5^vi which can be identified from a difference between 
(Ji Ji)^''*((5) and (J2</2)^''*(Q) ^^ the same Q^, for instance. 

Figured] shows the numerical results for B({Q), B2{Q)/a^ and C(i{Q)/a^ at the smallest 
quark mass (rra^ = 0.015) as a function of {aQY fo^' both vector and axial channels. In the 
momentum region {aQY < 2.3 only the B({Q) term gives sizable contribution, while the 
others are an order of magnitude smaller even without the suppression due to {aQY- Their 
dependence on {aQY is rather mild, so that it seems reasonable to fit these functions as a 
polynomial of (aQ)"^. We use a third-order polynomial to model these functions. This is 
used to subtract the artifacts at the momentum points for which the above procedure is not 
applicable, e.g. below the lowest (aQ)^ ~ 0.776. 

We notice that the difference between J = V and J = A is consistent with zero within 
statistical errors. This indicates that these lattice artifacts are strongly constrained by 
the exact chiral symmetry of the overlap fermion, and the effect of the finite quark mass 
is negligible. It also suggests that such short distance quantities are insensitive to the 



spontaneous chiral symmetry breaking, as it should be. This property is essential in the 
calculation of the difference II\MQ) — nV(Q), which is related to the electromagnetic mass 

nn 

difference of pions [18|, 119| . 

C. Perturbative calculation of the lattice artifacts 

Since the lattice artifacts are most significant in the high (aQ)^ region, perturbative 
analysis of the discretization effects is expected to give a reasonable estimate. We calculate 
the vacuum polarization functions in the lattice perturbation theory at one-loop level, which 
means that only the zeroth order of a^ is included. We then extract the terms corresponding 
to I\^^^'\Q), Bl,iQ) and C^,{Q). 

We calculate the vacuum polarization diagram in which two (axial-)vector currents ([3]) 
and (jlj) are inserted. The renormalization factor Z is set equal to 1 at this order. In the 
momentum space, the two-point function is written as 

(VM-iQ) - £|^Tr[(l-^Z.o(Al)^oW7. 

(8) 



^ Do{K - Q) ] So{K - Q)j, 



2mQ 
where the fermion propagator So{K) is given by 



So{K) = J- 
2mo 

with 



uj{K) + b{K) 



+ 1 



(9) 



LU 



[K) = Y.'in\K,) + h{K)\ (10) 

h{K) = ^(l-cos(Ag)-mo (11) 

M 

for the overlap fermion and Dq{K)^^ = Sq{K). We set a = 1 in this subsection. In the 
perturbative calculation, mo may be set equal to 1. At the perturbative level, the vector 
and axial- vector current correlators are equivalent in the massless limit, because of the exact 
chiral symmetry of the overlap fermion. 

After performing the numerical integral in (^ we extract Bi2(Q), Cu(Q) and UyiQ) 
in ([5]) through the same numerical procedure as we used in the non-perturbative extrac- 
tion. To be explicit, we take representative values of {aQY between 0.4 and 2.3 and 
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FIG. 2: Momentum dependence of B({Q), C('^{Q)/a? and B2{Q)/a? calculated in perturbation 
theory. 

consider two different momentum configurations aQ^^^ and aQ^'^\ The results for Bi{Q), 
B2(Q)/a^, and Cii{Q)/a'^ are shown in Figure [2l As we found in the non-perturbative 
calculation, the (aQ)"^ dependence is rather mild and we may precisely model these func- 
tions by quadratic functions: B[^{Q) = 0.06930(59) -0.00332(85)(aQ)2 + 0.00009(27)(ag)^ 
B^^iQ) = 0.0025(22) + 0.0023(30) (aQ)^ - 0.0009(9) (aQ)^ and C[^^iQ) = -0.00507(14) + 
0.00227(20) (aQ)2 - 0.00046(6) (aQ)^. The fit curves are shown in Figure [1 

The same curves are also plotted in Figure [1] by dashed lines. These perturbative results 
show reasonable agreement with the lattice data. It indicates that the lattice artifacts are 
indeed well described by the perturbation theory. 



D. Results for the vacuum polarization functions 

Lattice results for the vacuum polarization function 11 j (Q) for J = \^ are shown 
in Figure [31 The vacuum polarization function can be extracted from off-diagonal // 7^ z/ 
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FIG. 3: riy {Q) from ofF-diagonal [li^ v and diagonal ^ = v correlators with (lower panel) and 
without (upper panel) the subtraction of i^j^jlQ) ^^'^ ^^^{Q). 



(triangles) and from diagonal fi = u (circles) components. Upper and lower panels show 
the data before and after the subtraction of B:^^{Q) and Cmn{Q)'^ (Q) terms. Namely, for 
the upper panel, II j (Q) is identified with the formula ([5]) but without the B^{Q) and 
^mniQ) terms. As discussed above, raw lattice data of the diagonal components receive 
large contamination from B({Q) while the artifact for the off-diagonal components is much 
smaller (below 0.01). 

After the non-perturbative subtraction of 5/2 (Q) ^'^^ C^iiiQ)^ '^^ observe that the off- 
diagonal and diagonal components give consistent results. It strongly indicates that the 
higher order lattice artifacts are unimportant. We average the diagonal and off-diagonal 
data in the following analysis. 
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III. OPERATOR PRODUCT EXPANSION 



V and A channels 



We now discuss the fit of tlie lattice data to tlie OPE expression of the form [2 



n(o+i) 



,2 ,.2^ , "^V^J.^2 ,.2n , r^J ir^2A'mm) 



m = c+co{Q^fi') + j^cM,f^n+cM 



OPE^^ ' '"^ '^ ' g2-™v^ .^ > ' -qqy^ , Q, 



,2,((a./vr)GG) 



n(o+i) 



+ccGm '' "q[ ' ■ (12) 

Instead of directly treating the Adler function, we analyze its indefinite integral 
(g2). The coefficient functions Co{Q\^^'), CHQ^fx^), C^^{Q^) and Cgg(Q') 
are analytically calculated in perturbation theory. The terms of order 1/Q^ and higher are 
not included. 

A constant c is divergent and thus scheme-dependent, while other terms are finite and 
well-defined. Although we need to specify the renormalization scheme, the scheme depen- 
dence should disappear as the higher order terms are included. The following formulae are 
consistently given in the MS scheme, so that the strong coupling constant as(/i) is defined 
in this conventional scheme. 

The leading term Cq{Q'^ , ^^) is known to 0{a1) in the massless limit [Tl, [Sj as 

Co(Q',/i') = T7-^<|-+41n^ + ^^^ y-16C(3)+41n^ 



167r2 [3 Q2 ^ 

2 r 



+ 



"^(/i^ 



TT 



41927 3701 , /1658 

Nf - 

216 324 -^ V 9 



11^;) c(3) + i55c(5) 



365 11 , /^ 8, \^,^, /ll l^r\^ ^J^\^ A^^ 



(13) 



where Nf denotes the number of fiavors, and the zeta function is numerically given as 
C(3) = 1.20205 ■ ■ ■ , C(5) = 1.03692 ■ ■ ■ . For a finite quark mass there is a contribution of 
0{im? /Q^) with running mass m = m{fi). This term is represented by C;^{Q'^, fi"^), which is 
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also calculated to C(a^) as 



«s(Ai^ 
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13 
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3 /i^ N ;i2 



(14) 



12-121„^ 

55 






47- 
167r2 



«s(At^ 



TT 



128 



32C(3) 



(15) 



We ignore terms of 0{m'^) and higher. 

The OPE corrections of the form {0^^'^)/Q^ start from the dimension-four operators mqq 
and {as/'K)GG. Their Wilson coefficients C^ JQ"^) and Cgg{Q^) ^'^^ known to 0{al) and to 
C(as), respectively, as [21 1 



cF(g') 



a^(/i2 



99 



1 a,(/x2^ 



+/- 



TT 



C^gg(Q' 
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1 
12 



TT L 24 TT 

3 7r 3V7rMv24 

2^. 2 X ,,2 



116 - 4A^/) + (66 - ANf) In ^ 



7 



36 



iV. 



11 

T 



6 



iVjlln 






■(4C(3)-3 + ln^)+0/4 



1- 



lla.(Q) 



(16) 
(17) 



18 TT 

Here we note that the "gluon condensate" {{as/7T)GG) is defined only through the perturba- 
tive expression like (fT2l) . Due to an operator mixing with the identity operator, the operator 



(as/vr)GG contains a quartic power divergence that cannot be unambiguous 



within perturbation theory, which is known as the renormalon ambiguity 22|]. Therefore, 



y subtracted 



the term {{as/'n')GG) in ( IT2l) only has a meaning of a parameter in OPE, that may depend 
on the order of the perturbative expansion, for instance. 

The quark condensate (qq) is, on the other hand, well-defined in the massless limit, since 
it does not mix with lower dimensional operators, provided that the chiral symmetry is 



12 



preserved on the lattice. Power divergence may appear at finite quark mass as ma~'^. In 
the OPE formula flT2|) . it thus leads to a functional dependence ■m'^a~'^/Q^. Since the quark 
mass in the lattice unit is small (0.015-0.050) and (aQ)"^ is of 0(1) in our lattice setup, 
this divergent contribution is tiny (~ 0.1-0.2%). In fact, we do not find any significant m^ 
dependence in the lattice data. We therefore neglect this m^ dependence in the numerical 
analysis. 

B. V — A channel 

In addition to the individual vector and axial-vector correlators, we consider the V — A 
vacuum polarization function. For the difference Ily^^iQ) = ^v \Q) ~ ^a \Q)i ^^^ 
lattice data are more precise than the individual IIj {Q), so that the 1/Q^ and 1/Q® terms 
are also necessary: 



n(o+i) 

^W-A 



,p, m = ici - ct)m^, + {ci - ct^) m^^ 



+ (aG{^) + bQ{^i)\n^ + CQmgj— + ^. (18) 

In the V — A combination the coefficients C^ — C^ and C^ — C^ start at 0(a,^). The 



coefficients a6(/i) and bQ{n) contain dimension six operators Og and Oi as [9|. 



m 



aM = 2'iT{asOs){fi) + lA{alOs)ifi) + 2{alOi){fi), (19) 

5 

bM = -(a^08)(/i) + ^(«20i)(/i), (20) 

and the definition of these operators is given by 

(Oi) = Yl ((9"«7MAV^.g,.)fe7MAVgg,) - {q^l,l,X^4q,Ml,l,X''4q,)) , (22) 



^l,a,^,J 



with generator matrices r^ and A" of flavor SU(2) and color SU(3) symmetries, respectively. 
The numerical coefficients in the definition of a^ and b^ correspond to those of the Naive 
Dimensional Regularization (NDR) of 75. 

Unlike the dimension-four quark condensate {mqq), (Og) and (Oi) remain finite in the 
massless limit, hence gives leading contribution. The term ce, which has a mass-dimension 
five, describes their dependence on the quark mass. The term a^/Q^ represents the contri- 
butions from dimension-eight operators. 
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IV. FITTING RESULTS 



A. Fit parameters 



In the fitting of tfie lattice data witli tlie functions (IT^ and flTSl) . we fix tlie scale /j, 
to 2 GeV. We use the value of the quark condensate obtained from a simulation in the 



e-regime using the same lattice formulation at slightly smal 



-[0.251(7)(11) GeV]= 



er lattice spacing, (gg)(2 GeV) 



24 



25| are slightly different from 



23l |. (The values quoted in [12, 
but consistent with this number. The precise value does not affect the fit much, since 
the contribution of the Cqg term is sub- dominant.) The quark mass is renormalized in 
the MS scheme using the non-perturbative matching factor ^^(2 GeV) = 0.838(17) [I4I as 
m{fi) = Zm{fi)fnq. The coupling constant asifJ^) is transformed to the scale of two-flavor 



QCD, A^, using the four-loop formula [26 1 



as(/i^ 



TT 



4 



1- 



/3i\nL 



-i^Sc"^— -)-!} 



^o'LnK*'"'^^^"''^''"^ 






(23) 



with 



Pi 
P2 



1 

4 
1 

42 
1 

43 

1 

44 



11 --Nf 
3 ^ 



102- 
2857 



38,, 

—Nf 
3 ^ 

5033 



2 
149753 
6 
/ 50065 



Nf 



325 ,,2 
54 ^ 



(24) 
(25) 
(26) 



3564C(3) 



1078361 6508 



I- 



162 



1?C(3),.JH- 



162 
1093 . 

m 

729 ^ 



27 



C(3) Nf 



(27) 



and L = Hf^yA^}') 



MS 



Then, the free parameters in the fit are the scheme-dependent constant c, the gluon 
condensate paemeter {{as/ 7^)00), and the QCD scale A^ for the fit of an average 



n 



MS 

y^^iQ) ^ nV"-^^'(g) +nr^^(g). For the difference U'^l'jiQ), A^ 



(0+1) , 



r{o+i). 



(0+1), 



,(2) 



obtained above is 



used as an input and the dimension-six condensates ag, b^ and Cg are free parameters. 
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FIG. 4: Fit range dependence of A^^^, {(as/7r)GG) and the constant term c. The maximum 
momentum squared {aQ)l^ax is 1.324. 



B. V + A channel 



The OPE analysis requires a window in Q^ where the systematic errors are under control. 
The upper hmit (aQ)max — 1-3238 is set by taking the points where different definitions of 
the lattice momentum, i.e. Q^ = (2/a) sin(7rn^/L^) and Q^ = {2/a)'Jin^/L^, give consistent 
results within one standard deviation. In the physical unit, this corresponds to 1.92 GeV. 
To determine (a<5)niin! '^^ investigate the dependence of the fit parameters on (aQ)min i^ 

I— 5-1 ('?\ 

Figure HI We observe that the results for ^jJq, {{as/7i)GG), and c are stable between 
(aQ)^;^ ^ 0.48 and 0.65, which correspond to the momentum scale 1.16-1.35 GeV. Above 
(aQ)^i^ ~ 0.65 the fit becomes unstable; the results are still consistent within one standard 
deviation. 

Figure [5] shows the lattice data for Hy,,,^ (Q) at each quark mass and corresponding 
fit curves. It is clear that the Q^ dependence of the lattice data is well reproduced by 
the analytic formula. The quark mass dependence of Hy^^ (Q) is, on the other hand, not 
substantial as expected from the fit function (TT2|) . Our fit with the known value of (qq) 
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FIG. 5: Ily I ^ (Q) as a function of {aQ)^. The lattice data at different quark masses are shown by 
open symbols. Fit curves for each quark mass and in the chiral limit are drawn. The full result in 
the chiral limit (dashed-dots curves are at the finite masses, and solid curve is in the chiral limit), 
as well as that without [asG"^) /Q^ term (dashed curve), are shown. 

reproduces the data well. In the chiral limit, flT^ is controlled by two parameters: ^-J-q and 
{{as/'iT)GG) (apart from the unphysical constant term c). The fit result in the chiral limit 
is drawn by a solid curve. The dashed curve drifting upwards towards low Q^ region shows 
the result when the contribution from the {{as/'n')GG) term is omitted by hand. It indicates 
that the Q^ dependence is mainly controlled by the perturbative piece while the dimension- 
four term gives a minor contribution, which becomes slightly more important in the low Q^ 
regime. Numerically, we obtain A^ = 0.234(9) GeV and {{aJn)GG) = -0.058(7) GeV^ 
from a global fit of the lattice data at four different quark masses. The fit range of (aQ)"^ is 
[0.65, 1.3238]. 

Figure [6] shows Ai^ extracted from the lattice data at each quark mass. The flat be- 
havior provides another evidence that the lattice data are consistent with the perturbative 
prediction (IT^ . 
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FIG. 6: A^-TTT from the data at each quark mass. 
MS ^ 



C. Systematic errors 



In this subsection we discuss on possible systematic errors in this determination of A 



(2^ 
MS' 



That includes an estimate of the discretization effects and that of the truncation of pertur- 
bative and operator product expansions. 

As indicated from the perturbative analysis presented in Section IllCt the discretization 
effects are estimated reasonably well using the perturbation theory. Here, we discuss on the 
one-loop results for Ilj{Q) on the lattice for our choice of the fermion action and the current 
operators. This aims at estimating the remaining systematic errors due to the discretization 
effects after explicitly subtracting the B;^{Q) and C;^n{Q) terms. 

We again calculate the same one-loop vacuum polarization diagram at representative 
values of {aQ)"^ between 0.1 and 2.0. After subtracting the B^2{Q) ^^d C^^{Q) terms 
determined perturbatively in Section III CI we numerically obtain the piece corresponding 
to Ilj{Q), which contains the physical logarithmic dependence — l/(47r)^ln((aQ^)) as well 
as the lattice artifacts. In the continuum theory (or the perturbative calculation with the 
dimensional regularization, to be specific) only this logarithmic term appears, hence we may 
identify the remaining terms as the lattice artifacts. They are parametrized by a polynomial 
of {aQf. 

The result of the one-loop calculation is shown in Figure [7] (upper panel). We fit the data 
with a function including the known logarithmic term plus a quadratic function of {aQY and 
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FIG. 7: One-loop calculation of Ilv{Q) (upper panel) at representative values of (aQ)'^ (circles) 
and a fit with the log-plus-polynomial form. The dashed curve shows the purely logarithmic 
contribution. The term that represents the lattice artifact AUviQ) is shown in the lower panel. 



obtain the numerical result n^^*P'^(Q2) = __1^ ln((aQ)2) + 0.03085(9) + 0.00952(30)(aQ)2 - 
0.00132(20)(aQ)l 

In order to estimate the impact of this size of the discretization effect, we add this term to 



the fit function ( fT2l) and repeat the whole analysis. The result is A^^ = 0.249(37) GeV and 



(2) 



MS 

{{as/Ti)GG) = +0.11(15) GeV^ We find that Ag^is not largely affected, while {{as/7i)GG) 
is very sensitive to the lattice artifact and in fact changes its sign. 

Other (Lorentz- violating) discretization effects due to B:^{Q) and C:^^{Q) are subtracted 
non-perturbatively so that the associated error should be negligible. With our preliminary 
calculation of the above mentioned conserved vector and axial- vector currents for the overlap 
fermion, we confirmed that the results are consistent with the calculation presented in this 
paper obtained with the non-conserved currents ([3]) and (j4j) up to the unphysical constant 
term c. This observation confirms that our procedure to subtract the B;^{Q) and G^^{Q) 
terms is working as expected. 

The truncation of the perturbative and operator product expansions is also a possible 
source of the systematic error. In order to estimate the size of the former, we repeat the 
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FIG. 8: Fit range dependence of aQ(fi), bQ(fi) and ag. The horizontal axis denotes the minimum 
momentum squared (aQ)^;^. 



analysis using the fit formulae truncated at a lower order (two-loop level), and find that the 

(2) 

change of ^-rrq is niuch less than one standard deviation. It indicates that the higher order 
effects are negligible. The error from the truncation of OPE is estimated by dropping the 
terms of 0{1/Q'^) from flT2|) . From fits with higher (a(5)mm (between 0.79 and 0.89) to avoid 
contamination from the 1/Q^ effects, we obtain Ai^ = 0.247(3) GeV. The deviation of A^^ 
is about the same size as that due to the discretization effect. 

The errors due to finite physical volume and the fixed topological charge in our simulation 



271] are unimportant for the short-distance quantities considered in this work. A simple order 



counting gives an error of order 1/{QL)^ < 0.4% or smaller. 

To quote the final result, we take the central value from the fit without the discretization 
effect 

A|^ = 0.234(9) (+^0') GeV, (28) 



where the first error is statistical and the second is systematic due to the discretization and 
truncation errors. The result is compatible with previous calculations of ag in two-flavor 

and 0.249(16) (25) GeV ^ . (The physical scale is 



QCD: A 



(2^ 
MS 



0.250(16)(16) GeV 



normalized with an input tq = 0.49 fm.) 
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FIG. 9: Q i^y_j^{Q) as a function of (aQ) . The lattice data at different quark masses are shown 
by open symbols. Fit curves for each quark mass and in the chiral limit are drawn. 



D. V — A channel 

For the fit of the V — A vacuum polarization I[y_J{Q), we also examine the fit range 
dependence. In Figure [8] the fit parameters ag, b^ and og are shown as a function of (aQ)^;^ 
while fixing (a(5)max ^^ ^^^ same value 1.3238. We attempt to fit with (filled symbols) and 
without (open symbols) the a^/Q^ term in order to investigate how stable the results are 
against the change of the order of the 1/Q^ expansion. We find that the fit with a^/Q^ is 
stable down to (a<5)min — 0.46, while the other could not be extended below (a(5)mm — 0.58. 
The difference between filled and open symbols is marginal for a^ (circles), but too large to 
make a reliable prediction for feg (squares). To quote the results we set (aQ)min = 0.586 for 



bothnS^+J(Q) andn 



^'^'\Q). 



V-A 

In Figure [9], we plot Q^YVy_2{Q) as a function of {aQY for four different values of the 
quark mass irtq. The quark mass dependence is clearly observed. The main contribution 
comes from a dimension-six term c^mq/Q^, while the dimension-four term (jnqq) /Q'^ is sub- 
dominant (~ 20%), as its coefficient starts at 0{as)- In the chiral limit, there is a small but 
non-zero value remaining in Q^n^_^ |ope(Q^) as shown by a dashed curve in the plot. This 
is due to the four-quark condensates a^ and h^. 
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The four-quark condensate a^ obtained from Ii.y_J^{Q) is 

a6(2 GeV) = -0.0038(3) (1^^) GeV^ (29) 

where the first error is statistical. The second error represents an uncertainty due to the 
truncation of the 1/Q^ expansion. The central value is taken from the fit with a^/Q^ in 
f fTSj) and the error reflects the shift when this term is discarded. The result agrees with the 
previous phenomenological estimates —(0.003 ~ 0.009) GeV^ (30|. The other condensate is 
less stable; we obtain 66(2 GeV) = +0.0017(7) GeV^ or -0.0008(2) GeV^ with or without 
the 0{1/Q^) term, respectively. 

V. CONCLUSION 

Many of the lattice calculations to date have analyzed the two-point correlation functions 
to extract physical quantities such as the hadron mass spectra and decay constants. Usually 
the exponential fall-off of the correlator at large Euclidean time separation is used to isolate 
the ground state contribution. In this way, however, many interesting pieces of information 
are lost. They are in the short and middle distance regime where the perturbative analysis is 
also applicable. We use the two-point current correlators calculated on the lattice to extract 
the strong coupling constant with the help of the continuum perturbation theory and the 
operator product expansion. The recent work by Allison et al. has exploited |5| the similar 
idea and applied it to the charmonium correlator to extract the charm quark mass and the 
strong coupling constant. 

With the exact chiral symmetry realized by the overlap fermion formulation, the analysis 
of the lattice data is simplified. For the case of the vacuum polarizations, the continuum 
form of OPE may be applied without suffering from additional operator mixings, such as 
the additive renormalization of the operator gg, which appears in the Wilson-type fermion 
formulations. We also obtain the four-quark condensates a^ and fcg, which are relevant to 
the analysis of kaon decays 9||. 

In principle, our analysis does not require lattice perturbation theory, which is too com- 
plicated to carry out to the loop orders available in the continuum theory. But the pertur- 
bative calculation is still useful to estimate the discretization effects, which is well-described 
by perturbation theory in the asymptotic free theories. 
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The result for the strong couphng constant is compatible with previous lattice calcula- 



tions. The size of statistical and systematic errors is also comparable with them. An o 



bvious 



extension of this work is the calculation in 2+1-fiavor QCD, which is underway [3l|]- We 
also study the improvement of the analysis by using the conserved current for the overlap 
fermion formulation. 
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